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<N ; Abstract 

> 

cn ' We re-examine the evolution of neutrino asymmetries in several four neu- 



trino models. The first case involves the direct creation of Ly^ by V(, <-> Vg 



Q\ . oscillations. In the second case, we consider the mass hierarchy m^^ ^ 



"T-i/ , m^^ , m^^ where v^ ^-> Vg oscillations generate a large L^^ and some of this 



i-C \ asymmetry is converted into L^^ by v^ <-> v^, oscillations. We estimate the im- 

- ' plications for BBN for a range of cosmologically interesting drn^ values. The 

present paper improves on previous published work by taking into account 
,^ . the finite repopulation rate and the time dependence of the distortions to the 

neutrino momentum distributions. The treatment of chemical decoupling is 

also improved. 



I. INTRODUCTION 

If light sterile neutrinos exist, then this will lead to important implications for early 
Universe cosmology. This is because ordinary-sterile neutrino oscillations generate large 
neutrino asymmetries for a large range of parameters ||l]-§]. This is a generic feature of 
ordinary-sterile neutrino oscillations. 

The implications of this phenomena are quite model dependent. Various implications 
of this phenomena have been discussed in a number of previous papers for a number of 
interesting models motivated by the existing neutrino anomalies |0-§|,^|-§. For example in 
Refs. 0,^ it has been shown that the maximal u^^ -^ Vg oscillation solution to the atmospheric 
neutrino anomaly is consistent with a stringent BBN bound of N^Pr^ ~ 3.6 (and may also 
be consistent with A^ij-/^ < 3 depending on the model PJ^). This consistency requires 

rriy^ ~ few eV (for \5m1^^^^J\ ^ 3 x 10~^eV^) thus placing the Vr in the interesting hot dark 
matter range. Of course this is also of great interest to short-baseline experimentalists! 

Of particular concern to this paper, is the 'low temperature' evolution of neutrino asym- 
metries which also has important implications for BBN. In Ref. [Q, we discussed the four 
neutrino model with m^^ ^ irtj^ ^my^^irty^. In this case a large L^^ asymmetry is generated 
by Vj. ^-^ Ug oscillations, some of which is transferred to L,^^ by z/,- <-> z/g oscillations. This 
has important implications for BBN since it allows NfP^ < 3, with N^P^ ^ 2.5 for a large 
range of parameters if L^^ > 0. Qualitatively similar results occur for other sterile neutrino 
models as has been shown in a number of recent papers P,H . One point of all this is that 
N^j^j^ < 3 is a serious possibility if light effectively sterile neutrinos exist. 

In the paper Ref. we considered the case of large \5m?\ where the Ly^ was created above 
1.5-2.0 MeV, so that its implications for BBN could be approximately discussed by using 
thermal neutrino distributions (i.e. the neutrino asymmetry was distributed with chemical 
potentials) which were approximately constant during the era when the neutron/proton 
ratio was changing (i.e. for T ~ 2 MeV). We also briefly estimated the effects for the direct 
production of L^^ by z/g -^ Vg oscillations, (this case is only possible if the other neutrinos 
are lighter or do not oscillate into the sterile neutrino). 

Recently, Ref. 0| has been criticised in Ref. where it is claimed that the time de- 
pendence of the neutrino asymmetry and finite repopulation rate (which was assumed to be 
instantaneous in Ref. [0] for temperatures above about 1.5 MeV) is of critical importance. 
Ref. also similarly criticises Ref. (which studied a quite different 4 neutrino model 
with approximately degenerate z/^, Vr states) but this is clearly unjustified because Ref. 
takes into account the finite repopulation rate using a Pauli-Boltzman approach (as well as 
the time dependence of the distortion in the neutrino distributions). In fact, Ref. [^j appears 
to follow the repopulation procedure advocated in Ref. |0 and re-examines the cases in Ref. 
using this repopulation procedure. In view of this, we have also decided to revisit the 
models considered in Ref. in this paper because we believe that the results of Ref. to 
be incorrect. We will compute the evolution of the number distributions taking into account 
the finite repopulation rate and time dependence of the asymmetry. As already emphasised 
above, such an approach was already used in the papers PJ^ discussing different models, so 
it is straightforward to apply it here. We will also improve on Ref. f^ by discussing more 
completely the effects of the two similar oscillation modes, z/^ ^-^ z^^, z/^ ^^ z/g. We also give 



a more accurate treatment of the kinetic decoupling region which suggests a shghtly lower 
kinetic decoupling temperature. 

II. BIG BANG NUCLEOSYNTHESIS 

The primordial deuterium to hydrogen (D/H) ratio can be used to give a sensitive 
determination of the baryon to photon ratio rj which, given the estimated primordial ^He 
mass fraction, can be used to infer the effective number of light neutrino flavours N^P/^ 
during the BBN epoch. This value can then be compared with the predictions for N^rj 
from various models of particle physics to find out which ones are compatible with standard 
BBN. For example, the minimal standard model predicts N^Pr^ = 3. At the present time, 
most estimates favour N^^^^ < 3.6 and some estimates favour A^ij-/^ < 3.0 |10|. Of course. 



even if a model of particle physics is shown to be incompatible with BBN, this does not 
necessarily mean that the model is incorrect, since for example, it is also possible that one 
of the standard assumptions of BBN may not be correct |]Tl] . 



For gauge models with effectively sterile neutrinos, one in general expects N^P^ ^ 3. In 
fact, N^P^ may be less than three or greater than three. The prediction for N^P^ depends 
on the oscillation parameters in a given model and also the sign of the asymmetry (which for 
various reasons cannot be predicted at the moment). One possible consequence of ordinary- 
sterile neutrino oscillations is the excitation of sterile neutrino states, which typically leads to 
an increase in the expansion rate of the universe and thereby also increases N^P^ . Another 
possible consequence of ordinary-sterile neutrino oscillations is the dynamical generation 
of an electron-neutrino asymmetry. This also has important implications for BBN, as it 
directly affects the reaction rates which determine the neutron to proton {n/p) ratio just 
before nucleosynthesis. If the electron neutrino asymmetry is positive then it will decrease 
N^P^ , while if it is negative then it will increase N^P^ . 

The neutron to nucleon ratio, X„(t), is related to the primordial Helium mass fraction, 
^p,byQ 

Yp = 2X„ (1) 

just before nucleosynthesis. The evolution of X„(t) is governed by the equation, 

fix 

_^ = _A(n -. p)X„ + \{p ^n){l- X„), (2) 

where the reaction rates are approximately 

\{n -^ p) = \{n + z/g — > p + e^) + \{n + e^ — > p + z/g) + A(n -^ p + e~ + u^,), 

\{p ^> n) = \{p + e~ — > ra + z/g) + \{p + T'e ^ n + e^) + \{p + e~ +Ve ^ n). (3) 

These reaction rates depend on the momentum distributions of the species involved. The 
2-body processes in Eq.(BD for determining n <->■ p are only important for temperatures above 



^For a review of helium synthesis, see for example Ref. |12| 



about 0.4 MeV. Below this temperature these weak interaction rates freeze out and neutron 
decay becomes the dominant factor affecting the n/p ratio. For example, an excess of v^ over 
I7e, due to the creation of a positive L^^ would change the rates for the processes in Eq.(|). 
The effect of this would be to reduce the n/p ratio, and hence reduce Yp. Neutron decay is 
not significantly altered by lepton asymmetries unless they are very large. It is quite well 
known that a small change in Yp due to the modification of Ve and V^ distributions does not 
impact significantly on the other light element abundances (see for example Ref. [T^). A 
small modification to the expansion rate, using the convenient unit N^Pr^ , primarily affects 
only Yp, with WM 



5Yp ^ 0.012 X 5N^/j^. (4) 

Since Appendix A of Ref. P], describes in detail how we compute the effect on Yp due to the 
modified Vf. and V^ distributions, we will not waste any more paper by discussing it further 
here. 



III. CASE 1: IMPLICATIONS FOR BBN OF ve ^ us OSCILLATION 

GENERATED L^^ 

In this section we will study the direct production of a large L^^ from Ue <-> z/g oscillations. 
We will ignore oscillations involving z/^ or z/^. This is only an approximately valid thing to 
do provided that either their masses are very small (so that the largest |(5m^| belongs to the 
t'e ^-^ '^s oscillations and the other oscillations have \6m'^\ much less than 1 eV'^) or that 
they do not mix with the z/g, Us (i-e. the z/g, Us decouple from the u^, Ur in the neutrino mass 
matrix) . 

Let us begin with some necessary preliminaries. Our notation/convention for ordinary- 
sterile neutrino two state mixing is as follows. The weak and sterile eigenstates Va (« = 
e,n,T) and Us are linear combinations of two mass eigenstates z/^ and z/^, 

Ua = cos Oasl^a + siu 9asl^b, Us = - siu Oasl'a + COS Oagl'b, (5) 

where O^s is the vacuum mixing angle. We define O^s so that cos 29^8 > and we adopt the 
convention that Srri^g = ml — m^. 

Recall that the a-type neutrino asymmetry is defined by 

In the above equation, n^ is the number density of photons, n^ = 2({3)T^/tt'^. 

Note that when we refer to "neutrinos" , sometimes we will mean neutrinos and/or an- 
tineutrinos. We hope the correct meaning will be clear from context. Also, if neutrinos are 
Majorana particles, then technically they are their own antiparticle. Thus, when we refer 
to "antineutrinos" we obviously mean the right-handed helicity state in this case. 

In Ref. [^ it was shown that ordinary-sterile neutrino oscillations generate large neutrino 
asymmetries for a wide range of parameters. This work built upon earlier work on ordinary- 



sterile neutrino oscillations in the early Universe [ITB]. As already discussed in detail in 



previous publications 0,|,0 the evolution of lepton number can be separated into three 
distinct phases. At high temperatures the oscillations are damped and evolve so that L^""^ — >■ 
(where L'-"-' = L^,^ + L^^ + L,^ + L^^ + r], and r] is related to the baryon asymmetry). In 
this region the resonance momentum for neutrino oscillations is approximately the same as 
anti-neutrino oscillations. If 5m^^ < then at a certain temperature, Tc, which is given 
roughly by 0, 

r,~i6( -^'";-;r"'°j 'MeV. (7) 

exponential growth of neutrino asymmetry occurs (which typically generates a neutrino 
asymmetry of order 10~^ at T ~ Tc, see figure 1 of Ref. 0] for some typical examples). Tak- 
ing for definiteness that the L,y^ is positive, the anti-neutrino oscillation resonance moves to 
very low values of p/T ~ 0.5 while the neutrino oscillation resonance moves to high values 
p/T ~ 10 (see Ref. [^ for a figure illustrating this). The subsequent evolution of neutrino 
asymmetries, which is dominated by adiabatic MSW transitions of the antineutrinos, follows 
an orderly 1/T^ behaviour until the resonance has passed through the entire distribution. 
The final asymmetry generated is typically in the range 0.23 ~ L,^^ ~ 0.37 [Q. Because 
the oscillations are dominated by adiabatic MSW behaviour it is possible to use a relatively 
simple and accurate formalism to describe the evolution of the system at the 'low temper- 
atures', T ~ Tc/2. In fact, we only need to know the values of the oscillation resonance 
momentum at T ^ Tc/2. Previous numerical work has already shown [^] that by T ^ Tc/2, 
neutrino asymmetry is generated such that 0.2 ~ p/T ~ 0.8 (the precise value depends 
on siii^ 6 as, Sm'^g). Furthermore the subsequent evolution is approximately insensitive to 
the initial value of p/T in this range (provided, of course, that negligible number of sterile 
neutrinos were produced at high temperature). 

In this section we will deal with the case of u^ ^^ Vg oscillations directly producing 
the Li,^ asymmetry. For the implications for BBN we are primarily interested in the "low 
temperature" evolution of the number distributions and lepton numbers in this case. Our 
analysis can be broken up into the following steps: 

1. We assume complete adiabatic MSW conversion of neutrinos at the MSW resonance. 

2. From this we can compute the evolution of lepton number asymmetries which not only 
dictates the momentum of the MSW resonances, but also the chemical potentials. 

3. Using these chemical potentials we can evaluate the equilibrium distributions from 
which we can estimate the actual distributions by a Pauli-Boltzman repopulation equa- 
tion. 

We now discuss each of these steps in detail. 

Consider, for the moment, two-flavour small angle (i.e. cos29cs — 1) ordinary-sterile 
^e "^ ^s neutrino oscillations. As discussed in detail in earlier papers the Vf. ^-> z/g neutrino 
oscillations only generate L^^ provided that ^ml^ < and this will be assumed in the 



forthcoming discussion^. We know from numerical integration of the exact quantum kinetic 
equations [^ that the adiabatic approximation is vahd provided that sin^ 2^es ~ few x 10~^°. 
Now, coherent small angle adiabatic MSW transitions completely convert z/g ^^ Vg at the 
resonance momentum of these states. The resonance momentum is given approximately by 
(see e.g. Ref. H), 



where Oq = A^/2C{'i)G p / t^"^ and 






L^''^ = 2L,^ + L,^. (9) 



In the above equation we have neglected the L,y^ asymmetry (as well as the baryon/electron 
asymmetry). This is because these asymmetries are unimportant in the low temperature 
region unless they happen to be large (i.e. greater than about 10~^). For adiabatic two- 
flavour neutrino oscillations in the early universe it is quite easy to see that the rate of 
change of lepton number is governed by the simple equation H], 






d{prcs/T) 



dT 



(10) 



where pres is the MSW resonance momentum, Eq.(P) and 



Xip) = -[N,M-N,M]- (11) 



n^ 



Also, Ny^{p) and N^^p) are the momentum distributions of the z/g and Vs states. In the 
above equations, the case L^^ > has been considered (so that the resonance occurs for 
antineutrinos). Equation (0) relates the rate of change of lepton number to the speed 
of the resonance momentum through the neutrino distribution. Reference [§] provides a 
detailed discussion of how this equation can be derived from the quantum kinetic equations 
for the case of adiabatic evolution with a narrow resonance width. As discussed in Ref. 0, 
Equation (|T0|) can be simplified using 

dT dT dL^^ dT ' ^ ^ 

from which it follows that 



dL,^ /X^(2^ -4/Xp, 



dT 1 _ fx^^Si2iIIl T^ + 



(13) 



dL^, ~^ L<^^ 



^ Note that |57Tig^| < m^ (where m,y^ is the mass of the state which is predominately fg). Recall 
that there is an experimental upper bound on rrii,^ which is a few eV if fg is a Dirac neutrino and 
about an eV if fg is a Major ana neutrino [16|. 
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where / = 1 for d{p^cs/T)/dt > (that is for d{p,-cs/T)/dT < 0) and / = — 1 for 
d{Pves/T)/dt < and we have dropped the momentum dependence of X in the above 
equation for notational clarity^. Eq.(|T3|) allows us to compute the evolution of L^^. As 
discussed earlier, it is valid from T ~ Tc/2 (with Pres/T ~ 0.5 at this point). For the more 
complicated multi-flavour case considered in section IV, coupled equations based on Eq.(|13|) 
will be used. We now must describe how we compute the evolution of N^^ . 
The MSW transitions effect the adiabatic conversion 



We) ^ Ws)- 

This means that as Pi sweeps through the V^ momentum distribution. 



(14) 



NvAPi) 



NvAPi). 
NvAPi)- 



(15) 



In our numerical work the continuous momentum distribution for each flavour is replaced 
by a finite number of 'cells' on a logarithmically spaced mesh. As the momentum Pi passes 
a cell, the number density in the cell is modified according to Eq. (P^ . Of course weak 
interactions will repopulate some of these cells as they thermalise the neutrino momentum 
distributions. The repopulation can also generate small chemical potentials for the other 
flavours (as will be discussed latter). We take repopulation into account with the rate 
equation for each flavour a = e, /x, r. 



d KM 

dtNo{p,T) 
dtNo{p,T) 



r„(p) 
r„(p) 



No{p,T) 



No{p,T) 
NvM 



Noip,T) Noip,T) 



(16) 



where Ta{p) is the total collision rate and is approximately given by 



Taip) = y.GlT'' 



P 



3.15T 



(17) 



with Hf. ~ 4.0,1/^^^ ^ 2.9 and Gp is the Fermi constant. Also, in Eq.(|l^), 
Nq{p,T), N'^'^[p,T, fiu^) are the equilibrium distributions with zero chemical potential and 
chemical potential fi^^ respectively^. Previous papers ||^J^ used a simple approximation 
whereby the transition out of chemical equilibrium occurred at the decoupling temperature 
'Pdec- Obviously this is not a sharp transition. Also, there will be small chemical potentials 
created by the other flavours as they create Fgi/g pairs to compensate for the loss of V^ states. 
In the appendix we discuss a more accurate (but more complicated) formalism to compute 



^ At T = Tc/2, / = 1 and it does not change sign during subsequent evolution 
^ Our convention for the sign of the chemical potential is N'^'^{p,T, fi^ "* — 



■«'' ~ 2^ l+ef.P+i^<^c)/T- 



the values of the chemical potentials for all of the flavours as a function of time (or equiv- 
alently temperature). The conclusion is that the simple treatment of chemical decoupling, 
discussed in previous papers is roughly valid, but noticable differences (though typically not 
greater than about 5N^Pf^ ~ 0.2) for BBN can occur. [Although there will not be much 
difference for the z/g <-> Vs case since the lepton number is generated so late, since experi- 
mentally |5mgg| ~ 10 eV"^ (see footnote 2)]. We will use the more complicated formalism 
discussed in the appendix to evaluate the chemical potentials for all of our numerical work 
in this paper. 

Using the above procedure the lepton asymmetry L^^, and the neutrino distributions, 
Ny^ (p, t), N-p^ (p, t) can be obtained. We can feed the N^^ (p, t), N^^ (p, t) distributions into a 
nucleosynthesis code (which we integrate concurrently) in order to compute the implications 
for BBN. It is useful to separate the total contribution to 5Yp into two contributions, 

5Yp = 6iYp + 62YP, (18) 

where ^lYp is the change due to the effect of the modified electron neutrino momentum 
distributions on the nuclear reaction rates, and 62YP is due to the change in the energy 
density (or equivalently the change in the expansion rate of the universe). While BBN 
is only sensitive (to a good approximation) to the total contribution, 6Yp, the separate 
parts will have quite different implications for the forthcoming precision measurements of 
the anisotropy of the cosmic microwave background. In particular it may be possible to 
estimate the expansion rate of the Universe at the time of photon decoupling [|l^] . 

The contribution 6iYp can be determined by numerically integrating the rate equations 
for the processes given in Eq.(^ using the modified electron neutrino momentum distribu- 
tions N^^ and N^^^ as discussed in Appendix A of Ref. 0]. The contribution 62YP can be 
computed from the momentum distributions of the ordinary and sterile neutrinos through 




a=l 



pdp-3], (19) 



where 



00 



po^j^ N,ip,T)pdp=^^T\ (20) 

is the energy density of a Weyl fermion at equilibrium with zero chemical potential. [Recall 
that Eq.(^) can be used to express 6Yp, 6iYp and S2YP in terms of effective neutrino number, 
6N^jr^^ , 6iN^j^^ and 62N^jr^^ , respectively.] To calculate 62YP, we numerically determine 
the momentum distributions at T = 0.5 MeV. Because of the approximate kinetic decoupling 
of neutrinos for temperatures below about 3 MeV, large contributions]^ to 62YP, should they 
exist, must have been generated earlier. A temperature of 0.5 MeV is therefore a safe place 
to evaluate the final S2YP. 



^By 'large contributions' we mean 82^^^ ~ 0.10 



Recall that there is an ambiguity concerning the sign of the L^^ lepton asymmetry. We 
have considered the Ly^ > case above for definiteness, but L,^^ < is equally likely a 
priori. Previous work Ref. |0] has shown that the sign is fixed in the region where the 'static 
approximation' is valid. This approximation assumes that the the asymmetry evolution 
is dominated by collisions and is sufficiently smooth Q. Importantly it is generally valid 
in the region T c::^ Tc where the neutrino asymmetry is initially generated provided that 
sin^ 26* ~ (9 (10~^) for Sm^ ~ —10 eV"^. Thus in this region of parameters the sign is fixed. 
For sin^2^ ~ (9(10~^) numerical integration of the quantum kinetic equations (including 
the momentum distribution of the neutrinos) reveals 0] that the sign oscillates [] for a short 
period at T ~ Tg. In the parameter region where the sign oscillates (and possibly in some of 
the parameter region where it doesn't oscillate) the sign may not be fixed when fiuctuations 



are considered [^ . It may be possible for different regions of space to have different signs 
of the lepton number, as first suggested in Ref. [|l| . Whether this happens or not is an open 
question at the moment and will depend on the size of the fiuctuations present. 

In any case, even if the sign is fixed, it cannot be predicted as it depends on the initial 
values of the neutrino asymmetries (as well as the baryon asymmetry). For the negative L^y^ 
case, the roles of particles and anti-particles are reversed. One consequence of this is that 
the signs of all the other asymmetries are also reversed. The quantity 6iYp will obviously be 
significantly affected by this ambiguity in sign, while 62YP will not be affected at all. This 
means that we have two possible values for the overall change in the effective number of 
neutrino flavours during BBN. The results of the numerical work is presented in Figs. 1-3. In 
Figure 1 we show the evolution of L^^ for three examples, 6mlg/eV^ = —0.1, —1, —10. We 
emphasise that the evolution is approximately independent of sin^ 29es so long as sin^ 29es 
satisfles 




few X 10-^° ~ sin^ 2^^. ~ few x 10"^ I -^-^ ) , (21) 



where the lower bound comes from adiabaticity while the upper bound comes from the 
requirement that u^ ^^ Ug oscillations do not populate the sterile states at high temperatures 



^ In fact, it has been shown in Ref. ||5| that this approximation is equivalent to the adiabatic 
Uniit of the quantum kinetic equations in the region where collisions dominate the evolution of the 
neutrino asymmetry. 



^ Note that the recent study in Ref. ||l^, which neglects the neutrino momentum distribution 
arrives at quite different results. They find that the region where the sign doesn't oscillate is 
much smaller, roughly, sin^ 20 ~ O (10~^) for 6m? ~ — 10 eV"^. Qualitatively, this is very easy to 
understand. This reason is that in the average momentum toy- model, all of the neutrinos enter 
the MSW resonance at the same time which significantly enhances the rate at which neutrino 
asymmetry is created at T = T^.. The rapid creation of neutrino asymmetry reduces the region 
where the oscillations are adiabatic Q. Also, the paper in Ref. |^0| similarly assumes that all of the 
neutrinos have the same momentum, but (perhaps not suprisingly) obtain quite different results 
to Ref. |1|. 
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before the neutrino asymmetry is initially generated, i.e. for T r^ Tc (see e.g. Ref. 0). As 
explained earlier, we start the low temperature evolution at T ~ T'c/2 with the value of 
Pres/T ~ 0.5. Of course the full evolution from T ^ T^ to T — ;> can be obtained from 
numerical integration of the quantum kinetic equations (see e.g. Ref. [Q). However, for the 
applications considered in this paper only the low temperature evolution is required, which 
is why we use the much simpler formalism based on Eq. (|T3|) . 

The implications for BBN are shown in Figure 2,3. Figure 2 treats the Ly^ > case, 
while Fig. 3 displays the L^^ < case. As these figures show, the effect of the z/g •«-»• Ug 
oscillations on BBN is very significant and depends sensitively on the sign of the asymmetry 
and on the Smig value. We emphasise that our equations contain approximations. The most 



important are that the repopulation is handled approximately via Eq.(16). It is obviously 



difficult to estimate the size of this uncertaintly without computing repopulation exactly. 
Nevertheless we expect that this theoretical uncertainty is typically less than 5N^Pf^ ~ 0.2. 

It is evident that 52Yp is close to zero for the range of 5m1g considered. This can be 
approximately understood by noting that the generation of sterile states occurs below the 
kinetic decoupling temperature for z/g's. This means that the z/g states which have converted 
into sterile states are not significantly repopulated. It is interesting to observe that 52Yp 
is computed to be slightly negative for values of 5m?^g ~ few eV^. This feature can be 
qualitatively explained as follows. As p/T moves through the initial part of the spectrum 
the MSW transitions deplete the I7g's before the Ve spectrum is significantly distorted by 
the Ly^ asymmetry. However, by the time p/T moves to the higher momentum part (i.e. 
p/T ~ 4), the neutrino asymmetry is so large that the number of I^e states are significantly 
reduced (c.f. with the Fermi-Dirac distribution with zero chemical potential). For this 
reason the oscillations deplete more low momentum neutrinos than high momentum ones. 
The concurrent thermalisation of the neutrino distribution evidently reduces the average 
p/T per neutrino, and hence the energy density divided by T^ can be slightly reduced. Of 
course the temperature of the V^ states would be expected to increase a little, which in our 
approximation is neglected (much of this temperature increase would be absorbed by the 
other fiavours which are still in approximate thermal equilibrium down to temperatures 1-2 
MeV). 

Note that in the case of very small \Sm'^g\/eV'^ ~ 10~^, the dominant effect of the neutrino 
asymmetry in the L^^ > case is the modification to the Pauli blocking of the neutron decay. 
Ordinarily neutrino asymmetries lead to a negligible effect for neutron decay however in this 
case the effect is small (but not completely negligible) because the neutrino asymmetry is 
so large (L^^ ~ 0.37 for 6mlJeV'^ lO'^). 

Finally note that our results are quite different to the results of Ref. 0. We do not know 
why this is the case. Unfortunately, Ref. gives few details of how they computed the 
evolution of the neutrino asymmetry. Thus it is difficult to know whether the difference lies 
in the details of the asymmetry evolution or in the BBN code. 

This concludes our study of the direct production of Ly^ from z/g <-^ Us oscillations. We 
now consider the alternative, but more complicated case of the indirect production of L^^ 
from a large L^,^ asymmetry. 



10 



IV. CASE 2: IMPLICATIONS FOR BBN IN THE FOUR NEUTRINO SCENARIO 

WITH AU^ » AU^ , M^^ , M^^ 

We now discuss the second scenario of Ref . , that is a four neutrino model with 

m^^ > m^^ , m^^ , m^^ (22) 

We will first deviate from Ref. by considering the simpler case where the muon neutrino 
is ignored. Strictly, this would only be possible if the mixing angle satisfies sin^ 2^^,- ~ 
few X 10"^''. This case is simpler because there are then only two important oscillation 
modes, z/^ ^^ Ug and Ur ^^ Ve- Moreover these two oscillation modes always have quite 
different resonance momentum, so that they can each be described by two flavour oscillations. 
Later (in part B) we will consider the alternative case where sin^ 26'^^ ~ few x 10~^° and 
Vj. <-i> v^ oscillations are also important. 

A. Decoupled muon neutrino 

If we ignore the muon neutrino then there are two oscillation modes with approximately 
the same |5m^|, which we denote as Sm^^^^^: 

Vr ^ T^s, v^ ^ z/e. (23) 

Note that ^rrii^ ~ m^^ given Eg. (|2^) . We will consider the parameter space region where 
the Sm"^ values of all the other oscillation modes are small enough so that they can be 
approximately neglected for temperatures T ~ 0.4 MeV. (this will be true if the 6171^ of 
these other oscillation modes are all much less than about leV"^). This last condition means 
that these modes will not significantly affect the neutron/proton ratio and hence cannot 
significantly affect BBN. 

In the following discussion we consider the case L^^ > for definiteness. This means 
that the 17^ ^^ Vg generate L^^ while Vj. ^^ V^. oscillations reprocess some of this asymmetry 
into Ly^. Note that in this scenario the sign of Ly^ is necessarily the same as the sign of Ly^. 

The evolution of this system can be described by a straightforward generalisation of the 
two- flavour case given in Eq. (p!3D . In this case there are two MSW resonances, V^ <-»■ Vs 
and Vr ^^ Ve- We denote the resonance momentum of these two oscillations by Pi and P2 
respectively. They are related to the neutrino asymmetries and temperature through the 
equations, 

T ~ aoTUi' T ~ aoTUs' ^ ' 

where Oq = Ay/2C{'i)G p / t^"^ andQ 



^ We neglect the Ly^ asymmetry (and small baryon/electron asymmetries) which is a valid thing 
to do provided that it is less than about 10~^. 
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Li = 2Lp^ + Lp^, L2 = Ly^ — Ly^. 



(25) 



The evolution of the lepton numbers can be obtained by a straightforward generahsation to 
Eqs.([Tg,|n|,[T2|,|T3D, 



dU 



dT 

It 

dL^ 
dT 



where 



X2 
0, 

T 



d{Pi/T) 



dT 
d{P2/T) 



Xo 



d{P2/T) 



dT 



dT 



X,^ — [N^^{P,)-N^XPi 



n^ 



X2^-[N^AP2)-N^XP2)]. 



n^ 



Expanding out Eq.(^) we find 



dL,. 



yi 



a + (3 



dL^^ 



where 



dT "^ dT 

dLy^ 5- dLy 



- 1 fY ^(-Pl/^) f Y ^(-P^/^) 1 ^ 2/iXiPi /2X2P2 



2/2 = 1 + /2X2 

a = fiXi 
P = fiXi- 
6 = -/2X2 



d{P2/T) 



d{Pi/T) 

dT 

d{Pi/T) 

d{P2/T) 



1 



+ /2X 



/2X2P2 
TL2 
d{P2/T) _ 
dT 



TL, 



TLo 



-4/1X1P1/T2-4/2X2P2/T2 



+ /2X 



9(P2/r) -/iXiPi , /2X2P2 



aL, 



TLi 



TLo 



4/2X2P2/r2 



p = -/2X2 



dT 

d{P2/T) /2X2P2 



9L„ 



PL. 



(26) 



(27) 



(28) 



(29) 



and fi = 1 for d{Pi/T)/dt > and /i = -1 for d{P,/T)/dt < {i = 1,2). Solving Eq.(|| 
we find, 

c^Li,^ 6/3 + y2a 



dT 2/22/1 - pP ' 



dL„^ 6 + p- 



dT 



2/2 
12 



(30) 



In order to integrate these equations we need to specify the values of L^^ (or equivalently 
Pi/T) at T = Tc/20. The high temperature evolution typically does note generate significant 
L^^, i.e. typically L^^ < L^^. So that we have Pi/T ~ 0.5, P2/T ~ 2Pi/T. 

The evolution of the number densities is treated in a similar fashion to the previous 
section. Specifically, the MSW transitions effect the adiabatic conversions 

Pr) ^ Ps), Pr) ^ Pe), (31) 

at p = Pi and p = P2 respectively. This means that as Pi, P2 sweeps through the 17,-, z/g 
momentum distributions, 

iV^^(Pi)^iV^^(Pi), 
iV^,(Pi)^iV^,(Pi), 

iV^^(P2)^iV^,(P2), 

N^^iP^) ^ N^^P^). (32) 

In our numerical work the continuous momentum distribution for each flavour is replaced 
by a finite number of 'cells' on a logarithmically spaced mesh. As the momentum Pi, P2 
passes a cell, the number density in the cell is modified according to Eq. (|3^) . Of course weak 
interactions will repopulate some of these cells as they thermalise the neutrino momentum 
distributions. We take re-population into account approximately with rate equations of the 
form Eq.(p!6|) (for both u^ = Ve and z/^, = i/^) and compute the chemical potentials via the 
procedure discussed in the appendix. 

As in the case of the previous section, the evolution is approximately independent of 
sin^26'rs, sin^ 26'Te so long as 

/ gT/2 \ \ 

few X 10"^° ~ sin^ 20„~ few x 10"^ ^ , 

few X 10"^° ~ sin^ 2^^^, (33) 

where the lower bounds comes from adiabaticity while the upper bound comes from the re- 
quirement that v-r *^ Us oscillations do not populate the sterile states at high temperatures 
before the neutrino asymmetry is initially generated. Of course we have implicitly assumed 
that the z/g <-> i/g, z/^ ^-> Ug oscillation modes do not significantly populate the sterile neu- 
trinos before BBN. This is a valid assumption for a large range of parmaters even if the i^s 
mixes with large mixing angles with z/^ and/or z/g. For example, if z/^ <-*> Ug oscillations are 
approximately maximal with |5?ti^^| ^ 3 x 10~^ eV'^ (as suggested by the atmospheric neu- 
trino anomaly) then the z/^ ^^ Vg oscillations can potentially populate the sterile neutrinos 
at the high temperature T ?« 8 MeV . However this doesn't happen if the z/^ is in the eV 
mass range (for a large range of sin^ 2^^^) [0)9- 



^ Of course we also need to specify the initial signs, /i, which we take as positive. Subsequent 
evolution does not change these signs. 
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In Figure 4 we plot the evolution of the neutrino asymmetries, L^^ , L^^ , for three exam- 
ples, ^^farge/^^^ = ^-^i 10) 1000. As these figures show, for low values of ^^farge/^^^ the 
transfer of L^^ from L^^ is not very efficient. This is expected because the transfer of asym- 
metry relies on repopulation to distribute the L^^ away from the momentum region where it 
is created (i.e. p ~ Pi ) to the momentum region where Vr <-> V^ oscillations are important 
(i.e. at p ~ P2). As Smf^.^^^ increases the temperature where the L^^ asymmetry is created 
increases which makes the transfer to L^^ more efficient because of the faster repopulation 
rate. 

The implications for BBN are shown in Figures 5,6. Figure 5 treats the L^^ > case, 
while Fig. 6 displays the L^^ < case. As before, our results have a theoretical uncertainty 
which is dominated by the approximate treatment of repopulation. This uncertainty is 
expected to be typically less than about 6N^fW ~ 0.2 (obviously we expect this uncertainty 
to be much smaller than this when 5iN^^j is small). In the L^^ > case there is a 
dip at around 5m^ ~ —15 eV'^. This can be qualitatively understood as follows. For 
6171'^ ~ —15 eV^ the repopulation rate is not so rapid. This has two obvious effects: First, 
the thermalisation of the Ur momentum distribution is not so rapid and this would lead 
to less efficient production of L,y^. This effect would lead to a decrease in L^^ (and hence 
decrease \6N^j-^'^\). Second, the effect of a slow repopulation rate on the z/g distribution 
would be expected to have the opposite effect. The reason is that the momentum at which 
Lu^ creation is most significant is in the high momentum tail. This is because P2/T ~ 2Pi/T 
and significant L^^ is not generated until L,^^ is sufficiently large, i.e. roughly Pi/T ~ 2. 
The distortion of the z/g distribution in the high momentum tail greatly enhances the effects 
for BBN because these effects depend quadratically on the neutrino momentum. Evidently 
our numerical work indicates that the second effect dominates over the first effect. 

Finally, we now consider the more interesting, but more complicated case with the muon 
neutrino included. 



When the muon neutrino is included (i.e. sin^ 29f^r ~ few x 10 ^^) there are three 



B. Including the muon neutrino. 

> 
oscillation modes with approximately the same |(5?7i^|, which we again denote as Stuy^^ : 

f^ ^ Us, Ur ^ t^t,, l^T ^ t^e- (34) 

All the other oscillation modes have much smaller \Sm'^\ values. As before, we will consider 
the parameter space region where the dm"^ values of all the other oscillation modes are small 
enough so that they can be approximately neglected for temperatures T ~ 0.4 MeV. This 
last condition means that these modes will not affect the neutron/proton ratio and hence 
cannot significantly affect BBN. 

In the following discussion we again consider the case L^,^ > for definiteness. This 
means that the Vr ^^ Vs generate Ly^ while the other two oscillation modes reprocess some 
of this asymmetry into L^^, L^^. In Ref. [|] this system was first considered in this context. 
There, it was assumed that L^^ = L^^. While this is a good approximation for large enough 
values of |(5m^|, it is not always valid (as we will show below). In the following we will not 
assume this and consider the effect of the three oscillations modes. 
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At this point one may legitimately worry about 3-flavour effects. This is because the 
resonance momentum of the Vr ^^ ^e and v^ <->■ i^^ oscillation modes are expected to be 
approximately equal. However it turns out that these oscillations tend to be dynamically 
driven apart as we will explain latter on. Thus, it turns out that it is actually reasonable to 
treat all three oscillation modes independently as 2-fiavour MSW transitions. In the earlier 
paper this issue was not fully discussed so our treatment here improves on Ref. 0]. 

In this system there are three MSW resonances, 17^ <-> z7s, 17^- 4-^ V^ and Vr ^^ Ve. We 
denote the resonance momentum of these three oscillations by Pi, P2 and P3 respectively. 
They are related to the neutrino asymmetries and temperature through the equations, 



""^large 



T aoT^Li 



(35) 



where z = 1, 2, 3, uq = AV2C{S)GF/7r^ and 



Li = 2L^^ + Lj, + L,y^ , L2 = Ly^ — Ly , L3 



L,. 



(36) 



Using a simpler procedure to the previous (sub) sections, we have 



dT 

dT 

dLu,, 
dT 



-Xi 
X3 

X2 



d{Pi/T) 



dT 
diPs/T) 



X,, 



d{P2/T) 



dT 



-X, 



d{P,/T) 



dT 



dT 
d{P2/T) 



dT 



(37) 



where 



Xi 



T 



n-. 



[Nj,APi)-NvXPi)], 



N^AP2)-N^iP2 



T 



X,^-[N^^{P,)-N^XP^)]. 



n. 



(38) 



Expanding out Eq.(^) we find 



dL,^ 



yi 



2/2- 



dT 
dT 



a + (3 



dT 



+ 7- 



dL. 



dT 



2/3- 



d_L^ 
dT 



dL^^ 



r] + 



dT 



(39) 



where 
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yi 






1 + 






dL„ 



2/2 = 1 + /3X3 
2/3 = 1 + /2X2 



5(A/T) 



9(P2/T) 



a = fiX 



dL, 

d{Pi/T) 



1 + 
1 + 



TL3 

/3^3-P3 
TL3 ' 

/2X2P2 



+ /2X, 



TL2 

9(P2/T) 



+ /3X 



9(P3/T) 



iX,P^/T' - 4/2X2P2/T2 - 4/3X3P3/T2, 

—fiXiPi f^XsP^ 



d{P,/T) ^ d{P,/T) 
P = Ji^i — 57 H /3A3 — — = 

., _ , ^ d{Pi/T) , „ ^ a(P2/T) 

^- -/3X3^^ = 4/3X3P3/^^ 



TLi 



TU 



-flXiPi /2-^2-P2 



TLi 



TL, 



P = -/3X 



dT 
diPs/T) hX^P, 



dL„ 



TL, 



-hX,^^^ = ^hX,P,lT\ 



= -/2X 



dT 

d{P2/T) /2X2P2 



dL,^ 



TU 



(40) 



and fi = 1 for d{P,/T)/dt > and /i = -1 for d{P,/T)/dt <0{i = 1,2,3). Solving Eq.i 
we find, 

o^^^e Sy-iiyiy-i - 70) + py-iio^ys + 777) 



dT 



2/22/3(2/12/3 - 70) - pPvi 

«2/3 + 7^ + /?2/3-''^""^ 



dT 



1 
2/3 



2/12/3 - 70 



(41) 



In order to integrate these equations we need to specify the values of L^^ (or equivalently 
Pi/T) at T = Tc/2. The high temperature evolution typically does not generate significant 
L^,,L^^ (i.e. typically L^„ L^^ < L^J. So that we have Pi/T ~ 0.5, P2/T ^ P3/T ~ 
2P1/T. Now, there is no reason why P2/T should exactly coincide with P3/T (although 
it will be approximately equal). In fact for these two oscillations, P2/T = P3/T is not 
dynamically stable as we shall now explain. The behaviour of oscillations such as these (in 
the context of a quite different model) has been studied in some detail in Ref . . Generically 
there are two possible outcomes. Either the evolution of lepton number is such that it drives 
P2/T -^ P3/T or the evolution of lepton numbers is such as to drive them apart. To figure 
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out whats going to happen imagine that P2/T is shghtly less than P^/T. In other words the 
Vt- <-> Ve oscillation resonance preceeds the 17^ ^^ 17^ oscillation resonance. This means that 
the Vt- <-^ z/g resonance will efficiently interchange u^- and 17^ states at the resonance. This 
will transfer some L^^ to L^^ and will thus speed up the resonance a bit since it is inversely 
proportional to the difference of L,y^ and L^^. The trailing z7^ ^-> 17^ resonance will be less 
effective in transferring L,y^ to L^, because at this resonance there will be approximately 
equal number of 17^ and Vr states thanks to the efforts of the 17^- <-*• !7g resonance. Thus 
the two resonances will slowly move apart until eventually they will be far enough apart 
so that the thermalisation due to the collisions will be rapid enough to thermalise the i^r 
spectrum such that the L^, is created at approximately the same rate as the L^^. For 
definiteness, we will assume as our initial condition that P3/T > P2/T. For our numerical 
work we will assume that P3/T = I.OIP2/T initially. This could be due to a slightly larger 
|(5m^g| > I^TTi^^l for example. We found very similar results for even smaller choices such as 
P3/T = I.OOIP2/T. 

It is straightforward to numerically integrate the evolution equations, Eg . (PD , with T = 
Tc/2 'initial conditions' as described above. We keep track of the number distributions 
of all the 4 flavours (using a completely analogous procedure to the previous cases). In 
Figure 7 we plot the evolution of the neutrino asymmetries, L^^, L^^, L^^, for three examples, 
"^Warge/^^^ = 0.1, 10, 1000. The implications for BBN are shown in Figures 8,9. Figure 8 
treats the L^^ > case, while Fig. 9 displays the L^^ < case. As these figures show, this 
case is very similar to the previous results were the muon neutrino was neglected. Of course 
if we had started with P2/T > P3/T then the L,^^ and L^,^ are approximately interchanged. 
In this case, the modification to N^Pf^ would be somewhat smaller, especially for lower 

values of (^^farge- 

Let us compare our results with the original work in Ref. ||^ and the more recent work 
of Ref. 0. In Ref. we made the approximation that the repopulation was instantaneous 
above about 1.5 MeV. We also assumed that L„^ = L^^ and derived evolution equations 
consistent with that assumption. We found that 6NP^'^ ^ —0.5 for L^^ > and 6NP^^ ^ 
0.4 for L^^ < for 

10 ~ 6ml^^JeV^ ~ 1000. (42) 

For the L^^ > case the results are in rough agreement for Smf^j. ~ 100 eV'^. The difference 
for larger ^mf^rge is due to the more accurate treatment of chemical decoupling which suggests 
a lower decoupling temperature. For the L^^ < case the effect is underestimated by about 
0.4 in 6NP^'^ in Ref. 0. This difference is partly due to a mistake in the numerical work 
of Ref. P) which we have recently discovered. 

In Ref. they consider the case of part A of this section, i.e. neglecting the muon 
neutrino. Their results seems do not seem to be consistent with ours, expecially for L,y^ > 0. 
We do not know the reason for this. 



V. CONCLUSION 

We have made a detailed study of several 'four neutrino scenarios' which can generate 
significant L^y^ asymmetry thereby affecting BBN (these scenarios were first discussed in this 
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context in Ref. 0|). In the first case we considered the direct production oi L,^^ from z/g ^-> Ug 
oscillations. Our results are shown in Figure 2,3. Clearly very large modifications to BBN 
are possible and depend sensitively on 6m'^ [but are approximately independent of sin^ 29 
so long as sin^ 26' is in the range, Eq.fplj)]. The results also depends critically on the sign 
of L^^. We also studied the indirect production of L^^ from L^^. Our results are given in 
figures 5,6 for the case where the muon neutrino can be neglected and figures 8,9 where the 
muon neutrino is included. These results are in rough agreement with our earlier conclusion 
0] that 6N^j^^ ~ —0.5 for the case of positive L^^ for the parameter range Eq.(^). Notice 
that the figures show a slightly larger effect for 6m'^ ~ 10 eV"^ where the slow repopulation 
rate becomes important (Ref. assumed that repopulation was instantaneous). Also the 
more accurate treatment of repopulation in the kinetic decoupling region suggests a larger 
5N^Pr^ for \5m^\/eV'^ ~ 100 than was found previously. 

We conclude by emphasising once more that the detailed predictions of models with light 
sterile neutrinos are quite model dependent. Quantitatively different results occur for four 
neutrino models with approximately degenerate i/^, v^ 0, as well as in six neutrino models 
with three light sterile neutrinos approximately maximally mixed with each of the ordinary 
neutrinos [§. It is a remarkable prospect that accurate determinations of the primordial 
element abundances, as well as forthcoming precision measurements of the anisotropy of the 
cosmic microwave background may one day help to distinguish between competing models 
of particle physics. 



APPENDIX A: REPOPULATION 

Consider for definiteness the case of V^ ^^ T's oscillations which generate a large L^,^. 
The Va <->■ T's oscillations depelete the V^ states at the MSW resonance. Elastic collisions 
will tend to thermalise the momentum distributions so that they can be approximately 
described by chemical potentials, while inelastic collisions will create and modify the chemical 
potentials. Let us denote the total elastic and inelastic collision rates by the notation, 
rf , r^ respectively. It happens that Ff ^ r^H ^^ i^ makes sense to describe the neutrino 
distributions in terms of chemical potentials and a common temperature throughout the 
chemical decoupling period (2 ~ T/MeV ~ 4). We emphasise that the actual momentum 
distribution is always computed from Eq. I^T^, the purpose of this present discussion is to 
work out the evolution of the chemical potentials which are needed in the right-hand side of 
Eq.(^. 

The chemical potentials are related to the lepton number by the equation: 

_ 1 /■°° x'^dx 1 r°° x'^dx 

"" ~ 4C(3) Jo 1 + e^+^" ~ 4C(3) Jo l + e^+^-' ^ '' 

where /zq, = fiu^/T and p^a = f'Uc/'^ ^^^ C(3) — 1.202 is the Riemann zeta function of 3. 
Expanding out the above equation. 



Numerically Pf /F^ ~ 6.3, Tfj^/Vi ^ ~ 8.0 



24C(3) 



TT (i^a - jj^a) - 6(yU„ - fi^) In 2 + (/i: 



■^4 



/i|^ 



(A2) 



This is an exact equation for fta = —fia, otherwise it holds to a good approximation provided 
that \fla,a\ ~ 1. If we turn off inelastic collisions for a moment, and assume that the 
oscillation generated L^^ is positive then the oscillations generate a large jla- The evolution 
of fia due to the generation of L^^ (i.e. due to oscillations) can easily be obtained from 
Eq. 



dfJ^a 



dT 



24C(3) 



TT^ — 12 jlaln2 + 3/i^ 



dL^^ diic 



dT ' dT 



0. 



(A3) 



Now, let us turn on the inelastic collisions and see what happens. There are six inelastic 



processes, which we list below together with their thermally averaged interaction rates |T^ 



process 


rate 


(1) v^V^ ^ Z/^Z7^ 


n = Fo 


(2) Vt-Vj. ^->- z/gFg 


ri = Fo 


(3) u^TJ^ ^^ e+e" 


r| = (8a;2 - 4x + l)Fo 


(4) l^^.'Ui, ^ l^eT^e 


ri = Fo 


(5) u^Uf, ^ e+e" 


Tl = (8x2 _ 4a, + i)^Q 


(6) UeVe <->■ e+e" 


r^ = (8x2 ^ 4^ ^ 1)^^ 



In the above table, x = sin^ ^^ is the weak mixing angle (sin^ 6*^ ^ 0.23), and Fq = 

^j^ n-y ~ 0.13G%T^. Consider first the first process listed in the above table. This process 
will change the number of Vr^Vr, Vfi^v^j, states such that 



dn,. 



dt 



driT,^ 



(1) 



dt 



dn„ 



(1) 



dt 



dnr 



(1) 



dt 



(A4) 



(1) 



where the subscript ' Lx' denotes the contribution to the rate of change due to the process 
(1) in the table. The rate can be expressed approximately as follows. 



d{n^Jno) 



dt 



(1) 



1 
no Jo 



oo roo 



'KMNu,{p') - N^MNvAp')] cri{p,p')dpdp' (A5) 



For convenience we have normalised with respect to the number of neutrinos in a Fermi- 
Dirac distribution with zero chemical potential, no = -^:;^ — • 
approximation that 



Let us further make the useful 



N.^ip) ~ e-l'-Noip), N^M - e-^^Noip), iV,±(p) = Noip), 



(A6) 



where Nq{p) is the Fermi-Dirac distribution with zero chemical potential. Note that the e+ 
distributions have zero chemical potential due to the very rapid collisions with the back- 
ground photons (see e.g.Ref |]12|)- Thus, with the above approximation, Eq. ( |A5|) can be 
expressed in the simple form. 
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dVu^ 



dt 



(1) 



{v<^^Vu^ - VurVu.) r{, 



(A7) 



where rj^^^ = e^°' , Tfp^ = e^" . Clearly, we also have that 



dVur 




drju^ 


dVu, 




dVu^ 


dt 


(1) 


dt 


(1) ^' 


(1) 


dt 


dVu. 
dt 


(1) 


_ dvv^ 
dt 


= 0. 

(1) 







(1) 



(A8) 



A similar set of equations can be obtained for the other five inelastic processes. Putting this 
altogether, we have 



drjur 



dt 



dr]^ 



drjjy. 
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dt 
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drju. 
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drju. 
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dt 



j^dr], 



repop 



i=l 



dt 



(i) 



dt 



diji^. 



repop 



dt 



Arfr/, 



repop 



i=l 



dt 



(i) 



Of course the total rate of change of fia, /i« is given by 



(A9) 



dHa dfia 



dt dt 

dt dt 



+ 



+ 



Hi LliQ/^ 



dt 
djli 



dt 



repop 



repop 



(AlO) 



The above equation can be used to approximately compute the set of chemical potentials. 
Then using Eqs.(|I^) the evolution of the set of number distributions can be obtained. 

In Figure 10 we give the evolution of jla,a for some illustrative examples. In the Figures 
we consider the model discussed in section IV (part A) with the parameter choice 5mf^^^^ = 
10 eV"^ for Figure 10a and Smf^^.^^ = 200 eV"^ for Figure 10b. 

We have also compared the above repopulation procedure with the simpler procedure of 
a fixed decoupling temperature, TJ^^ where fi^^ = —fiva fo^ T > TJ^^ and fi^^ frozen (in the 
case L^^ > 0) for T < TJ^^. We get rough agreement provided that Tj^^ ~ 3 MeV. 
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Figure Captions 

Figure 1: Low temperature evolution of L^^/h (h = T^/T^) due to z/g ^^ ^s oscillations, 
for the parameter choices, 6m'^ = —0.1 eV'^ (solid line), 6m'^ = — 1 eV^ (dashed line) and 
6m'^ = —10 eV'^ (dashed-dotted line). 

Figure 2: 5N^j^^ versus |(5mg^|. The dashed line is the contribution (5iiV^j^ due to the 
effects of the L^^ asymmetry while the dashed-dotted line is the contribution 52N^^^ due 
to the change in the expansion rate. The solid line is the total contribution 6N^j-^^ = 
SiN^/J^ + 62N^ff^. This figure considers the case L^^ > 0. 

Figure 3: Same as Figure 2 except L^^ < is considered. 

Figure 4: Low temperature evolution of L^^/h, L^^/h {h = Tj^/T^) for the model of case 2 
(part A). Figs. 4a,4b,4c correspond to the parameter choices (^mfarge — O-^ ^V^, ^"^wge — 
10 eV"^ and Smf^^^^ = 1000 eV"^ respectively. 

Figure 5: SN^j^^'^ versus ^^farge ^^ ^^e model of case 2 (part A). The dashed line is the 
contribution (5iA^fj-j^ due to the effects of the L^^ asymmetry while the dashed-dotted line 
is the contribution 52N^Pf^ due to the change in the expansion rate. The solid line is the 
total contribution 5N^/j^ = ^liVfj-/^ + 62N^jj^ . This figure considers the case L^^ > 0. 

Figure 6: Same as Figure 5 except L,y^ < is considered. 

Figure 7: Low temperature evolution of L^^/h, L^^/h, L^^/h (h = Tj^/T^) for the model 
of case 2 (part B). Figs. 7a,7b,7c correspond to the parameter choices Smf^ = 0.1 eV"^, 
^"^fargc = 10 eV^ and Snii^j.^^ = 1000 eV"^ respectively. 

Figure 8: 5NfPf^ versus Srn^^ for the model of case 2 (part B). The dashed line is the 
contribution 5iN^M^ due to the effects of the Lj^^ asymmetry while the dashed-dotted line 
is the contribution 52N^M^ due to the change in the expansion rate. The solid line is the 
total contribution 5N^Pf^ = 6iN^P^ + 62N^P^ . This figure considers the case L„^ > 0. 

Figure 9: Same as Figure 8 except L^^ < is considered. 

Figure 10: Evolution of fla,a for the model of case 2 (part A). Fig 10a, 10b corresponds 
to the parameter choices Smf^j.^^ = 10 eV"^ and Smf^j.^^ = 200 eV'^ respectively. In the 
figures the thin solid, dashed, dashed-dotted, dotted and thick solid lines correspond to 
fte, fj'fj., P'T, P^e, P^T respectively. Note that ft-p: = ft^ in this case. 



22 



to 

CO 




T/MeV 



0.20 







5m'/eV' 







5m'/eV' 



C5 



IW 


; 


' ' 1 




; 




- 




Figure 4a 


- 


i 






L 


^Nv 








T 


\v 






10' 


: 


\ 




: 


10-^ 


r 


\ 




- 






"" - -. \ 








L 


~^-.^ \ 




- 




e 


s \ 




■ 












10-^ 


r 


\ \ 


. 


- 




- 


\ 


\ 


- 






^ 


\ 






- 


\ 
\ 
\ 


\ 


- 




- 


\ 


\ 


- 


in-4 




\ 
V 1 


\ 





10' 



10" 



10^ 



T/MeV 



■-J 




T/MeV 



00 




T/MeV 



0.60 

g^BBN 



CO 




5m'/eV' 



1.20 

g^BBN 



eff 



CO 
o 



1.00 



0.80 



0.60- 



0.40- 



0.20- 



-0.20 




10' , , 10' 

5m'/eV' 



CO 



IW 


; 


' ' 1 




; 




- 




Figure 7a 


- 


i 


^"^ 




L 


^\^ 








T 


\^^ 






10' 


: 


\ 




: 


10-^ 


r 


\ 




- 




L 


"~'~-^^.^ \ 




- 




e 


s \ 
\ \ 

\ \ 






10-^ 


r 


\ \ 

\ 


V 


- 




L 


\ 
- - ^ \ 


\ 


- 




^i 




\ 




in-4 




. .^^^ \. 


\ 





10' 



10" 



10^ 



T/MeV 



CO 



IW 


; 


1 ' ' 


: 


i 


- 




Figure 7b 




L 






1 






10' 


r 




- 




L 




■ 




e 




■ 






"""""■ ^s. 




10-' 




"" s. N. 


- 




: L 










\ 




- 


^ X 


\ 




. 


^^ \ 


\ 




- 




\ 


10-^ 


r 


\ \ 

\ \ 


\ : 




- 


\ \ 
\ \ 

\ V 

\ \ 

\ V 

\ \ 


V 


in-4 




1 . \ 


_1 1 1 1 1 1 1 



10' 



10" 



10^ 



T/MeV 



10' 



10' r 



CO 
w 



10"' r 



10" r 



10-' 



; 


; 


" 


Figure 7c 


L 


"■~\^ 


X 


N. - 


L 


\ 


e 


>v 


;__ 


\^ 


r L 




^^ 


** 

\ 




\ 


i 


\ 




\ 


■ 


♦ 




♦ 




* 




♦ 




♦ 




\ 


~ 


i 




* 


- 


% 




\ 


_ 


\ 




\ 




\ ' 




\ 




\ 











10' 



10" 



10' 



T/MeV 



g^BBN 



CO 




6mVeV' 



1.00 



erf 



CO 



0.80 



0.60- 



0.40 



0.20 




5m'/eV' 



00 

C5 




T/MeV 



la/T 



CO 



l.UU 


Figure 10b 


0.80 


_____^ 


0.60 


\ 


0.40 


\ 


0.20 


\ 




;,.,,,,,,,,,,,,...,:.^,..,:..,:..::,__^^^ ^----__ 



0.20 


r-. -" ' ' 



10" 



T/MeV 



10' 



